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1. Introduction 

The low-energy behavior of QCD, the spontaneous breaking of chiral symmetry, including 
the explicit breaking by the quark masses, is described by chiral perturbation theory (XVT). In the 
lattice regularization of QCD lattice artifacts can contribute to the breaking of chiral symmetry, 
for Wilson fermions, or its partial breaking, in the case of staggered fermions. These effects can 
be included in the ZPT approach, leading to new, lattice discretization dependent low-energy con- 
stants. We consider Wilson fermions in this contribution, for which the effective theory, Wilson 
ZPT (WZPT) was introduced and worked out in [p. The new terms in the chiral Lagrangian affect 
the low-lying spectrum of the (Hermitian) Wilson-Dirac operator [||]. For a recent review with 
additional references, see Ref. |^. Here, we test and verify the predictions for the distribution of 
the low-lying eigenvalues with lattice QCD simulations [Q, ||] and show that they can be used to 
obtain the new low-energy constants introduced in WZPT. We also demonstrate the effect of clover 
improving the Wilson-Dirac operator. 



2. The WZPT and Wilson RMT framework 

We will be concerned with the e -regime of WZPT where the zero momentum modes dominate 
- the system size is such that m^cL -C 1. In addition we adopt the power counting with m ^ a^. 
Hence, dropping the kinetic part of the chiral Lagrangian, we consider 

1 1 
^ = --mZTr [U + U^) - -zITr [U - U^) + a^Y . (2. 1) 

The second term, representing a 1//^ temi, is introduced for later convenience. Y describes the 
lattice artifacts [|l|] 

r = W^i:v{u'^ + U^'^) +W(,\j:r:{U + U^)]^ +Wi\j:v{U -U^)f . (2.2) 

At large Nc, the two-trace terms are suppressed. 
The finite size scaling considered is such that 

m = miy , z = zLV and dj = a^WjV for 7 = 6,7,8 

are held fixed. Here Z is the condensate and V the volume. 

This leading order in WZPT can equivalently be described by a chiral random matrix theory 
(RMT). For Wilson fermions, including the one-trace term with low-energy constant W^, the Dirac 
operator is represented in Wilson RMT (WRMT) as 



SA iW 



(2.3) 



with W a random (n + v) xn complex matrix, and A and B random Hermitian matrices of size 
{n + v) X {n + v) and nxn, respectively. As usual in the RMT context, we consider a fixed index 
V. We use a chiral basis with 75 = diag (1,...,1,— 1,...,— 1).A and B represent the chiral symmetry 
breaking term corresponding to the Wilson term in the Wilson-Dirac operator. 
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Ens 




ro/a 


a [fm] 


size 


L[fm] 


|e| = 0, l,2cfgs 


A 


2.635 


5.37 


0.093 


164 


1.5 


1279, 2257, 1530 


B 


2.635 


5.37 


0.093 


204 


1.9 


401, 682, 644 


C 


2.79 


6.70 


0.075 


204 


1.5 


1207, 2130, 1448 



Table 1: The ensembles used. The scale is set by ro = 0.5 fm. The r^/a values come from interpolation 
formulae in [p^. Q is the topological charge (see text). 



The two-trace terms can be incorporated in WRMT via two Gaussian integrations 

1 



Z\/ / ^ ^ A. ^ ^ \ 
(m,z;a6,fl'7,fl'8j 



1 T 

^6 _ 



dy(,dyie. 



Z''(in-y6,z-y7;0,0,a8) . 



Here, 



Z^(m,£;0,0,a8) =/^/t/det^[/e-^^(W5=w,=o) 
is the fixed-index partition function with the one-trace 0'{a^) term included. 



(2.4) 



(2.5) 



3. Index of the Wilson-Dirac operator 

As indicated above, RMT predictions apply to gauge field sectors with a fixed index, or, in the 
continuum, fixed topological charge. For the Wilson-Dirac operator, the index can be defined by 



V = £'sign((^j75j^)) 



(3.1) 



with \k) the ^'th eigenstate of the Wilson-Dirac operator, D-w. Only eigenvectors with real eigenval- 
ues contiibute, and the / indicates that only the real eigenvalues in the branch near zero, with eigen- 
values < Vcut, are kept. Introducing the Hermitian Wilson-Dirac operator D5{mo) = Y5{D]v + mo) 
and using 

D5(mo)|t/A) =0 ^ Dw\v) = -mo\¥) (3.2) 

the index can equivalently be obtained from the zero crossings of the spectral flow of D5{mo) up 
to nicut = —fcut 10]. It corresponds to the index of an overlap operator [||] with kernel D=,{mcut)- 
Because of the dependence on the choice of r^,,,, the index of the Wilson-Dirac operator is not 
unique. 



4. The numerical simulations 

For our numerical tests, in the quenched case, we generated three ensembles using the Iwasaki 
gauge action [^, which suppresses dislocations and gives a fairly unique index v or topological 
charge Q. The ensembles are characterized in Table |[ 

The topological charge listed in Table |] was obtained after six steps of HYP smearing [|lT]] with 
an improved lattice FF operator [12|. On the configurations with \Q\ < 1, as well as the \Q\ =2 



configurations of ensemble A, we also did the much more expensive computation of the index 
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Figure 1: Distribution of the real eigenvalues of the Wilson-Dirac operator for the v = 1 configurations of 
the two L = 1.5 fm ensembles A and C (left) and comparison of the distribution with and without clover 
improvement at a = 0.093 fm (ensemble A, right). 



from the spectral flow. We first applied one HYP smearing before constructing the Wilson-Dirac 
operator. The topological charge and the index agreed on most configurations, with the agreement 
improving at smaller lattice spacing and becoming worse for the larger volume ensemble C, for 
which it was about 97%. 

The crossing points in the spectral flow are the real eigenvalues, whose distribution is shown 
in Fig. |l| (left). The dashed vertical lines are estimates of (minus) the critical mass. Some real 
eigenvalues are smaller, on so-called "exceptional" configurations. 

For ensemble A, we also computed the spectral flow with clover improving the Wilson-Dirac 
operator, again after one HYP smearing. The clover coefficient was set to 1, which is expected to 
be close to the nonperturbative value after the HYP smearing [pj|]. The resulting distribution of the 
real eigenvalues is compared to the unimproved case in Fig. |l] (right). The improvement is quite 
dramatic, besides the expected reduced shift away from zero, the distribution is much narrower 
and more symmetric. The width is determined by the ff{a^) terms in Eq. (2.2), so with clover 
improvement the coefficients are much smaller. 



5. Wilson eigenvalue distributions and WRMT 



We next computed the lowest 20, in magnitude, eigenvalues of the Hermitian Wilson-Dirac 
operator D5(nio) with bare mass amo = —0.216 for ensembles A and B to compare to eigenvalues 
distributions obtained from WRMT [^. 

Without clover improvement, we considered only contributions from the two-trace term in 
Eq. (2.2) [Q]. We used the v = histogrammed eigenvalue distributions of ensemble A to determine 
the WRMT parameters m and d = and the eigenvalue rescaling factor LV. Using the same 
parameters we then get a prediction for the v = 1 distribution that can be compared to the numerical 
data (see Fig. ^ top). 

Ensemble B differs from ensemble A only in the volume. Using volume scaling, mg = 
'Ma (Vb /Va ) and ds = d^ ^/Vb/Va, we obtain predictions for the distributions for ensemble B (see 
Fig. |2| bottom). As can be seen, the WRMT predictions work well. 
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Figure 2: Comparison of the histogrammed eigenvalue distributions with WRMT. The v = distribution 
of ensemble A (top left) was used to obtain the parameters. The predictions for ensemble B (bottom) used 
volume scaling of the parameters. 



For ensemble C, at the smaller lattice spacing, we computed the eigenvalues with two different 
bare masses aniQ = —0.178 and —0.184. We used the histogrammed v = distribution with bare 
mass amo = —0.184 to the determine the WRMT parameters, and used "mass scaling". Am = 
AmoLV for predictions for the distributions with the other bare mass amo = —0.178, as shown in 
Fig. |3[ Again, the WRMT predictions work well. 



6. Clover improved eigenvalue distributions and WRMT 

We have already seen from the distribution of the real eigenvalues in Fig. |l| (right) that clover 
improvement not only, as expected, decreases the additive mass renormalization (the real eigen- 
value peak is much closer to zero) but also the size of the ff{a^) low-energy constants considerably 
(the distribution becomes much narrower). Here we consider the effects on the distribution of the 
20 lowest, in magnitude, eigenvalues of the Hermitian Wilson-Dirac operator D=,{rnQ) with clover 
improvement for ensemble A using a bare mass airiQ = —0.03. The comparison with WRMT is 
shown in the first three panels of Fig. ^ 

For \dj \ ^ 1 the lattice effects affect, to leading order, only the index peak of the topological 
modes. These are the lowest eigenvalues with almost chiral eigenvectors which correspond to the 
real eigenvalues shifted by the bare mass. As can be seen in Fig. ^ (bottom right) the distributions 
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Figure 3: Comparison of the histogrammed eigenvalue distributions with WRMT for ensemble C. The v = 
distribution with bare mass amo = —0.184 (top left) was used to obtain the parameters. The predictions for 
bare mass amo = —0.178 (bottom) used "mass scaling". Am = AmoLV. 



match almost perfectly. The eigenvalue density, for | v | > 0, on the opposite side of the index peak 
is almost continuum like and allows determination of m and LV. We use the |v| = 1 eigenvalue 
distribution for this. The aj are then obtained from their effect on the index peak. We use the fact 
that the low-energy constant have fixed signs W(, < 0, Wj < and Wg > [|2|, 14, |l^] and that the 
distribution depends only on the combination jWel + [VVyl [Q] allowing to take Wj = 0. We find that 
with either ag ^ ot ^ we can reproduce the histogrammed | V | = and | V | = 1 distributions 
in Fig. ^ (top) equally well. But only with d^^O can we reproduce the | v| = 2 distribution, too, as 
shown in Fig. |^ (bottom left). 

We can explain the drastically different effect of and Wg on the analytic prediction for 
|v| = 2 by noting that the Wg-term, in WRMT, coiTcsponds to a Gaussian fluctuating mass, see 
Eq. (2.4). The 5-function index peak of the continuum theory is therefore smeared into a Gaussian 
peak with an amplitude that increases with |v|. Wg, therefore, does not introduce a repulsion 
between eigenvalues. On the contrary, the Wg-term of WZPT is included in the representation of 
the Dirac operator, Eq. ( [2.3[ ), of WRMT, and hence induces an eigenvalue repulsion, as can be seen 
from the red curve in Fig. || (bottom left). This repulsion is seen for all sectors with |v| > 1 
It is thus useful to include eigenvalues from configurations with | v | > 1 for a determination of the 
low-energy constants from fits to eigenvalue distributions. 

We finally note that, with clover improvement the nonvanishing \d^\ is about a factor 3-4 
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Figure 4: Comparison of the histogrammed clover-improved eigenvalue distributions with WRMT for en- 
semble A (top and bottom left). The v = 1 distribution (top right) was used to obtain the WRMT parameters. 
The red curves are the WRMT predictions with ag 7^ 0, the blue curves those with a^ ^ 0. The bottom right 
plot shows a comparison of the distribution of the real eigenvalues with the first two positive eigenvalues of 
Dj{mQ) shifted by the bare mass for the V = 2 configurations. 



smaller than the nonvanishing \a^\ without the improvement, both after one HYP smearing, illus- 
trating again the quite dramatic effect of clover improvement on the ff{a^) low-energy constants. 

7. Conclusion 

We have presented numerical simulations, in the quenched case, of the low-lying eigenvalues 
of the Hermitian Wilson-Dirac operator, both with and without clover improvement to compare to 
predictions from £-regime Wilson XPT or, equivalently, Wilson RMT. We used the Iwasaki gauge 
action which suppresses dislocations and leads to a fairly unique index or topological charge. This 
is helpful, since the analytical predictions are made for sectors of fixed index. We found that our 
eigenvalue distributions agree well with the analytical predictions, and verified scaling with volume 
and (bare) mass. 

We have also looked at the distribution of the real eigenvalues of the Wilson-Dirac operator, 
obtained from the spectral flow. We found a dramatic decrease of both the additive mass renormal- 
ization (the real eigenvalue peak is closer to zero) and the ff{a^) low-energy constants (the width 
of the distribution becomes narrower and more symmetric) with the clover improvement. 
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Fits to the distribution of the low-lying eigenvalues of the Hermitian Wilson-Dirac operator al- 
low determination of the low-energy constants of QCD including those that parameterize the tff{a'^) 
lattice effects. However, distributions on configurations with | v | > 1 are needed to disentangle the 
effects of from those of and W-j when all \dj \ are small. 
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